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Abstract 

We investigate a simple model to explain inflation, neutrino masses and dark matter simultane- 
ously. This is based on the so-called radiative seesaw model proposed by Ma in order to explain 
neutrino masses and dark matter by introducing a Z2-odd isospin doublet scalar field and Z2-odd 
right-handed neutrinos. We study the possibility that the Higgs boson as well as neutral compo- 
nents of the Z2-odd scalar doublet field can satisfy conditions from slow-roll inflation and vacuum 
stability up to the inflation scale. We find that a part of parameter regions where these scalar 
fields can play a role of an inflaton is compatible with the current data from neutrino experiments 
and those of the dark matter abundance as well as the direct search results. A phenomenological 
consequence of this scenario results in a specific mass spectrum of scalar bosons, which can be 
tested at the LHC, the International Linear Collider and the Compact Linear Collider. 
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I. INTRODUCTION 



The new particle with the mass of 126 GeV which has been found at the LHC [l, 2] is 
showing various properties that the Higgs boson must have. It is likely that the particle 
is the Higgs boson. If this is the case, the Standard Model (SM) of elementary particles 
is confirmed its correctness not only in the gauge interaction sector but also in the sector 
of electroweak symmetry breaking. By the discovery of the Higgs boson, all the particle 
contents in the SM are completed. This means that we are standing on the new stage 
to search for new physics beyond the SM. There are several empirical reasons why we 
consider the new physics. Phenomena such as neutrino oscillation 3-81, existence of dark 

n un 

matter [9| and baryon asymmetry of the Universe I9 M111 cannot be explained in the SM. 
Cosmic inflation at the very early era of the Universe [ijj], which is a promising candidate 
to solve cosmological problems such as the horizon problem and the flatness problem, also 
requires the additional scalar boson, the inflaton. 

The determination of the Higgs boson mass at the LHC opens the door to directly explore 
the physics at very high scales. Assuming the SM with one Higgs doublet, the vacuum 
stability argument indicates that the model can be well defined only below the energy scale 
where the running coupling of the Higgs self-coupling becomes zero. For the Higgs boson 
mass to be 126 GeV with the top quark mass to be 173.1 GeV and the coupling for the 
strong force to be a s = 0.1184, the critical energy scale is estimated to be around 10 10 GeV 
by the NNLO calculation, although the uncertainty due to the values of the top quark mass 
and a s is not small 3]. The vacuum seems to be metastable when we assume that the 
model holds up to the Planck scale. This kind of analysis gives a strong constraint on the 
scenario of the Higgs inflation 14| where the Higgs boson works as an inflaton, because the 



inflation occurs at the energy scale where the vacuum stability is not guaranteed in the SM. 
Recently, a viable model for the Higgs inflation has been proposed, in which the Higgs sector 
is extended including an additional scalar doublet field |l5 |. 

In order to generate tiny masses of neutrinos, various kinds of models have been proposed. 
The simplest scenario is so called the seesaw mechanism, where the tiny neutrino masses 
are generated at the tree level by introducing very heavy particles, such as right-handed 



neutrinos 



161 ] , a complex triplet scalar field 17J , or a complex triplet fermion field 



18j. The 



radiative seesaw scenario is an alternative way to explain tiny neutrino masses, where they 
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TABLE I: Particle contents and their quantum charges. 



are radiatively induced at the one loop level or at the three loop level by introducing Z 2 -odd 
scalar fields and Z 2 -odd right-handed neutrinos 19l-l2l|. An interesting characteristic feature 
in these radiative seesaw models is that dark matter candidates automatically enter into the 
model because of the Z 2 parity. 

In this Letter, we discuss a simple model to explain inflation, neutrino masses and dark 
matter simultaneously, which is based on the simplest radiative seesaw model 20] . Both the 
Higgs boson and neutral components of the Z 2 -odd scalar doublet can satisfy conditions for 
slow-roll inflation 22] and vacuum stability up to the inflation scale. We find that a part of 



the parameter region where these scalar fields can play a role of the inflaton is compatible 
with the current data from neutrino experiments and those of the dark matter abundance as 
well as the direct search results 23(. A phenomenological consequence of scenario results in 
a specific mass spectrum of scalar fields, which can be tested at the LHC, the International 
Linear Collider (ILC) [2J| and the Compact Linear Collider (CLIC) [25]. 



II. LAGRANGIAN 



We consider the model, which is invariant under the unbroken discrete Z 2 symmetry, 
with the Z 2 -odd scalar doublet field <£> 2 and right-handed neutrino v R to the SM with the 

n 

SM Higgs doublet field $i [20|. Quantum charges of particles in the model are shown in 
Table [H Dirac Yukawa couplings of neutrinos are forbidden by the Z 2 symmetry. The Yukawa 
interaction for leptons is given by 

C-Yuka W a = YJ^^R + Yj7 L ^\v R + h.C, (1) 
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where the superscript c denotes the charge conjugation. The scalar potential is given by [15] 



V = ^ + (6l$i| 2 + 6l<£ 2 | 2 ) J R 

+/i 2 |$i| 2 + ^|$ 2 | 2 + ^A^il 4 + ^A 2 |$ 2 | 4 + A 3 |$i| 2 |$ 2 | 2 
+A 4 ($ t 1 $ 2 )($ t 2 $ 1 ) + [iA 5 (($ t 1 $ 2 ) 2 + h.c.)}, (2) 

where Mp is the Planck scale (Mp ~ 10 19 GeV), and is the Ricci scalar. 

We assume that /i 2 <0 and /i 2 > 0. $i obtains the vacuum expectation value (VEV) 
v (= a/— 2/U 2 / Ai ~ 246GeV), while $ 2 cannot get the VEV because of the unbroken Z 2 
symmetry. The lightest Z 2 -odd particle is stabilized by the Z 2 parity, and it can act as 
the dark matter as long as it is electrically neutral. The quartic coupling constants should 
satisfy the following constraints on the unbounded-from-below conditions at the tree level; 

A x > 0, A 2 > 0, A 3 + A 4 + A 5 + ^\~ 2 > 0. (3) 

Three Nambu-Goldstone bosons in the Higgs doublet field $i are absorbed by the Z and W 
bosons by the Higgs mechanism. 

Mass eigenstates of the scalar bosons are the SM-like Z 2 -even Higgs scalar boson (h), the 
Z 2 -odd CP-even scalar boson (if), the Z 2 -odd CP-odd scalar boson (A) and Z 2 -odd charged 



scalar bosons [H^). Masses of these scalar bosons are given by [20| 

m H = A + + A 4 + A 5 )f 2 , 

m A = A + + A 4 - A 5 )v 2 , 



m 2 H± = i4+ ^X 3 v 2 . (4) 



III. CONSTRAINT ON THE MODEL FROM INFLATION AND DARK MATTER 
A. Inflation 

We consider the Higgs inflation scenario [ljl 26 1 in our model defined in the previous 
section. The scalar potential is given in the Einstein frame by 

v „ A: + A 2 r 4 + 2(A 3 + A 4 )r 2 + 2A 5 r 2 cos(2g) / _ 2 ^\ 2 , . 

E ~ 8(6r 2 + 6) 2 V 1 " 6 J' (5) 



where 0, r and 6 are denned as 



1 v^L ' ^ Le* 



, 6^ 



2^ 1 + ^ + 



M 2 P 



h2 



(6) 



with taking a large field limit £Ji\/M%> + £ 2 hl/Mp > 1. 



For stabilizing r as a finite value, we need to impose following conditions 15]; 



A 2 6 - (A 3 + A 4 )6 > 0, 
A16 - (A3 + A 4 )6 > 0, 
AiA 2 - (A 3 + A 4 ) 2 > 0. 



(7) 



Parameters in the scalar potential should satisfy the constraint from the power spectrum 9|, 



« / 2(A 1 + q2A 2 -2q(A3 + A 4 )) „ 4 
s2W A 1 A 2 -(A 3 + A 4 )^ " 5X10 ' 



(8) 



^ «A 2 - (A 3 + A 4 ) < 12 

C 2 Xx + a 2 A 2 - 2a(A 3 + A 4 ) ~ ' 1 ; 

where a is given as a = £i/£ 2 . When the scalar potential satisfies the conditions in Eqs. 
OZJ)-©, the model could realize the inflation. 

B. Dark Matter 



We assume that the CP-odd boson A is the lightest Z 2 odd particle. (By changing the 
sign of the coupling constant A5, the similar discussion can be applied with the CP-even 
boson H to be the lightest.) When A5 is very small such as O(10~ 7 ), A is difficult to act 
as the dark matter because the scattering process AN — > HN opens, where iV is a nucleon. 
The cross section is too large to be consistent with the current direct search results for 
dark matter 27], |28] . In Ref . 15] , the authors claim that both the Higgs boson and Z 2 -odd 



neutral scalar bosons can work as the inflatons when the dark matter (H or A) has the mass 
of 600 GeV if A 5 < 10 -7 . However, as recently discussed in Ref. |28(, the bound from direct 



search results are getting stronger, and such a dark matter is not allowed anymore in this 
model without a fine tuning among the scalar self-coupling constants. Therefore, we here 
consider a 6-digit fine tuning, 



a\ 2 - (A 3 + A 4 ) ~ 10" 



(10) 



at the inflation scale 1 . In this case, A5 can be taken to be a sizable value. The scatter- 
ing process AN — > AN then comes mainly from the diagram of the SM-like Higgs boson 
mediation. The cross section is given by 29] 



a(AN AN) 



\ 2 

A hAA 



111 



N 



-f 2 

2 J N' 



(11) 



where Xhaa = A3 + A4 — A5, /V = )^ m /y . fr 7 — f m v /ty; and tun is the mass of nucleon, where 



fr u + frd = 0.056, fxs = 30] and Jtg — 0.944 [31] . The mass itia should be approximately 
a half of irih [32J in order for the dark matter to be consistent with the abundance from the 
WMAP experiment [9] and the upper bound on the scattering cross section for AN — > AN 
from the XENON100 experiment 23j. The coupling constant Xhaa should satisfy 



Xhaa < 0.036, 



(12) 



at the low energy scale for consistency to satisfy the data from the XENON100 experiment. 



C. Tiny Neutrino Masses 



In this model, tiny neutrino masses are generated by the one loop diagram in Fig. [TJ 20]. 



The neutrino mass {m v )ij are given by 



[m 



\ 5 v 2 {Y V )1{Y U 



Z_> 167r 2 
k 



M k R 



m 



H 



In 



111 



H 



11 1 



H 



{M k nY {M k R y 



In 



m t 
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(13) 



where M\ is the Majorana mass of v\. The flavor structure of (m u )ij is given by 
(Yj,)!? (Y v )j / Mft. The neutrino mixing matrix is explained by neutrino Yukawa cou- 
pling constants (Y u ) k . The magnitude of tiny neutrino masses can be explained when 



1 When we take an one order of magnitude fine tuning for a\2 — (A3 + A4) ~ 0, it would satisfy the condition 
from inflation and dark matter for m A ~ 1 TeV. However, we here do not consider these parameter sets 
because it would be difficult to test at the LHC and the ILC. 
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FIG. 1: The Feynman diagram for tiny neutrino masses. 



{Y v yl (Y^j/Mr ~ O(10~ n ) GeV -1 because A5 and masses of scalar bosons, m H and m A , 
are constrained from the conditions of inflation and dark matter. Our model then can be 
consistent with current experimental data for neutrinos {3-8]. 



D. Running of Scalar Coupling Constants 



In the SM, the energy scale cannot reach to the inflation scale because the quartic coupling 
constant of the Higgs boson is inconsistent with the unbounded-from-below condition at 
10 10 GeV scale when m t = 173.1 GeV and a s = 0.1184 [l^j]. On the other hand, if we 
consider extended Higgs sectors such as the two Higgs doublet model, the vacuum stability 
condition on the quartic coupling constant for the SM-like Higgs boson can be relaxed due 
to the effect of the additional quartic coupling constants [3^. Therefore, these models can 
be stable up to the inflation scale. We calculate these coupling constants by using the 
renormalization group equations with the following beta functions 34j ; 
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FIG. 2: Running of the scalar coupling constants. Red (solid), blue (dashed), brown (dot-dashed), 
green (doted) and black (long-dashed) curves show Ai, A2, A3, — A4 and A5, respectively. 
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16tt 2 
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6A!A 3 + 2AiA 4 + 6A 2 A 3 + 2A 2 A 4 + 4A 2 + 2A 2 + 2A 2 

2 



+ -/ + - A g* - \g 2 g' 2 - 3A 3 (3s 2 + g 2 ) + QX,yt 



(17) 



(18) 



/3(A 4 ) 
/3(A 5 ) 



16tt 2 
1 

16tt 2 



2A 4 (A X + A 2 + 4A 3 + 2A 4 ) + 8A 2 + 3g 2 g' 2 - 3X A (3g 2 + g' 2 ) + 6X 4 y 2 ] , (19) 



2A 5 (Ax + A 2 + 4A 3 + 6A 4 ) - 3A 5 (3^ 2 + g' 2 ) + 6A s2 / 2 



(20) 



We here impose the conditions of triviality 



Ai < 27T, 



(21) 
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Ai 


A 2 


A 3 


A 4 


A 5 


10 2 GeV 


0.26 


0.34 


0.51 


-0.49 


0.0047 


10 17 GeV 


1.6 


6.3 


6.3 


-3.1 


0.0061 



TABLE II: An example for the parameter set which satisfies constraints from the inflation and the 
dark matter at the scales of O(10 2 ) GeV and O(10 17 ) GeV. 



and vacuum stability (the unbounded-below-condition) up to the inflation scale. In Fig 121 
running of the scalar coupling constants are shown between the electroweak scale and the 
inflation scale. The vacuum instability due to Ai is avoided by the effect of the Higgs self- 
coupling constants with Z 2 -odd scalar bosons |33| . In Table [Til we show an example for the 
values of the scalar coupling constants at the scales of O(10 2 ) GeV and O(10 17 ) GeV, which 
satisfy the conditions of the inflation and the dark matter, where O(10 17 ) GeV denotes the 

nn 

inflation scale for £ x ~ £ 2 = O(10 4 ) [14j, [15( in our model 2 . 
E. Mass Spectrum 

Let us evaluate the mass spectrum of the model under the constraint from inflation, the 
neutrino data and the dark matter data as well as the vacuum stability condition. In our 
model, there are nine parameters in the scalar sector; i.e., £i, £ 2) ^2> A3, A 4 and 

Ab. 

First of all, as the numerical inputs, we take v = 246 GeV and = 126 GeV. Second, we 
use the conditions to explain the thermal fluctuation; i.e., the equation where ~ is replaced 
by = in Eq. (jSJ). We also impose the condition where ~ is replaced by = in Eq (fTUl) . 
Third, the allowed region for the mass of the dark matter A is determined from the current 
invisible decay results at the LHC and from the constraint from the direct search results 
for the dark matter. The former gives the lower bound and the latter does the upper 
bound; i.e., 63 GeV< m A < 66 GeV. We here take m A = 65 GeV as the reference value. 



2 When fx ~ £ 2 = O(10 4 ), unitarity is broken at M P /£i [35|. Then, we should introduce new particle at 



the unitarity breaking scale to save unitarity 36|. However, we do not consider the effect of this particle 
on the running of A coupling constants because this effect affect only above O(10 15 ) GeV. The effect is 
expected to be smaller than the effect of the SM particles. 
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Further numerical input comes from the perturbativity of A2 up to the inflation scale; i.e., 
--^(/Uinf) = 27r, where /z; n f is the inflation scale 10 17 GeV. The parameter set in Table HT1 can 
be consistent with these numerical inputs and the constraints are given in Eqs. (J3J), (J7|), (j2J), 
QT2"]) and OH]). The mass spectrum of the scalar bosons is determined as 

m h ~ 126 GeV, 
m H ± ~ 140 GeV, 
m H — 67.1 GeV, 

rru ~ 65.0 GeV. (22) 

The mass spectrum is not largely changed even if m A is varied with in its allowed region. 
Consequently, in our scenario, the following relation for the mass spectrum is obtained; 

m A < m H^ 69 GeV ) < m h(~ 126 GeV ) < m H±(~ 140 GeV )- ( 23 ) 

The upper bounds on m H and vnn± are obtained in order to satisfy the conditions from 
Eqs. ()3]), (J7j) and (ED)- Therefore, we can test the model by using the mass spectrum at 
collider experiments. 



IV. PHENOMENOLOGY 



Masses of Z 2 -odd scalar bosons have been constrained by the LEP experiment. In our 
scenario, rn H± should be around 140 GeV, which is above the lower bound given by the 



LEP experiment 
than m 



37 



23, 



38j . From the Z boson width measurement, m H + m A should be larger 



391 ] . In addition, there is a bound on HA production at the LEP. However, 



when rriv f — m A < 8 GeV, masses of neutral Z 2 -odd scalar bosons are not constrained by the 



LEP 



37 



39]. The contributions to the electroweak parameters [40|] from additional scalar 



bosons loops are given by 4l|, |42] 



AS = -— [F' A (m H ±,m H ±) -F' A (m H ,m A )] 



AT 



V2G } 



16n 2 aEM 



[-F A (m A ,m H ±) - F A (m H ,m H± ) + F A (m H ,m A )] 



(24) 
(25) 
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where 



Fa(x,v) = F A (y,x 

F' A (x,y) 



x 2 + y 2 



x 2 y 2 , x 2 
In — 

ar — y z y z 



F' A (y,x) = -- 



4 x 2 In x 2 — y 2 In 
3 x 2 — y 2 



x 2 + y 2 

,2A2 



(26) 
(27) 



(3r — 

In all of our parameters, it is consistent with current electroweak precision data with 90% 
Confidence Level (C.L.) [42]. 



The detectability of H, A and H ± at the LHC has been studied in Ref . [43| , where the 
center-of-mass energy -y/i = 14 TeV and the integrated luminosity L = 300 fb^ 1 are assumed. 
They conclude that it could be difficult to test pp — > AH + / H H + / H + H processes because 
cross sections of the background processes are very large. The process of pp — > AH could 
be tested with about the 3cr C.L. when inn and uia are 50 GeV and 80 GeV, respectively. 
However, it would be difficult to test pp —> AH in our model. In our parameter set, and 
rriA are about 65 GeV and 67 GeV. In this case, after imposing the basic cuts [43(, event 
number of pp — > AH is negligibly small. 

We now discuss signals of H, A and if ± at the ILC with = 500 GeV. In the following, 
we use Calchep 2.5.6 for numerical evaluation 44J. First, we focus on the H pair production 
process e + e" — > Z*(Y) ~^ H + H~ — > W +( ^W~^> AA —> jjivAA, where j denotes a hadron 
jet 45]. The final state of this process is a charged lepton and two jets with a missing 
momentum. The energy of the two-jet system Ejj satisfies the following equation because 
of the kinematical reason; 



m* H ± - m\ 



< E ii < 



m H± — m 



^s-2^s/A-m 2 H± 



Ejj is evaluated by using our parameter set as 



17 GeV < Ejj < 180 GeV. 



(28) 



(29) 



— >■ 



The distribution of Ejj of the cross section for e + e — > Z*(^y*) — > H + H 
W+(*)W~(*) AA — >• jjiuAA is shown in Fig. |3j The important background processes against 
e+e" ->• Z*(Y) ->■ H+H W + ^W~^AA jjivAA are e+e" W + W~ jjlv and 
e + e~ — > Z{j)Z — > jjii with a missing £ event. In these processes, the missing invariant 
mass is zero. These backgrounds could be well reduced by imposing an appropriate kine- 
matic cuts. We expect that m H ± and can be measured by using the endpoints of Ejj at 
the ILC after the background reduction. 
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1 




Second, we focus on HA production e + e~ — > Z* — > HA — > AAZ* — > AAjj at the 
ILC The final state is two jets with a missing momentum. Ejj is evaluated by using our 
parameter set as 

0.28 GeV < Ejj < 15 GeV. (30) 

The distribution of Ejj for e + e~ ->■ Z* -)■ ifA ->■ AAZ* ->■ AAjj is shown in Fig. U 
We note that the lower endpoint is nearly zero in this case because the mass difference 
between H and A is small in our parameter set. The main background processes against 
e + e~ — > Z* — > HA — > AAZ* — > AAjj are e + e~ — > ^Z — > ^jj with missing photon event 
and e + e~ — > ZZ — > jjvv. In the e + e~ — > ^Z — > ^jj process, Ejj distributes around 
i/s/2 because the initial energy is shared by the outgoing photon and the Z boson. In the 
e + e" — > Z(j)Z — > jjvv processes, the missing invariant mass distributes around the Z boson 
mass. These backgrounds could also be reduced by imposing kinematic cuts. Therefore, we 
can measure m# and ttla by observing the endpoints of Ejj at the ILC. 

Finally, we discuss prediction on the diphoton decay of the Higgs boson h. BR(h — > 77) 



in the model, which the SM with Z 2 -odd scalar doublet, has been studied in Ref. 



46] . The 
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FIG. 4: The distribution of Ejj for the cross section for e + e 



HA -> AAZ* 



AAjj. 



deviation in our model from the SM is given by 
BR(h 77) 



N c Q}A 1/2 (r t ) + A^tw) + M^A (t h ±] 



(31) 



BR(h SM -> 77) |iV c Q 2 A 1/2 (r t ) + ^(r^)! 2 

where N c and Qf are the color and electromagnetic charges of the top quark, respectively. 
A 1 / 2 (x), Ai(x) and A (x) denote 

A 1/2 (x) = 2{x + (x-l)f(x)}x- 2 , 

Ax(x) = -{2z 2 + 32; + 3(2x-l)/(x)}ar 2 , 

A (x) = -{x-f(x)}x- 2 , (32) 
where t x and f(x) are given by 



2m T 



In 



arcsin 2 (y / x) 

l+y/l-l/a 



27T 



X < 1 



X > 1 



When we use our parameter set in Eq. ( )22|) . the ratio is calculated as 

-> 77) 



BR(h SM -> 77) 



0.92. 



(33) 
(34) 



(35) 
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In our model, BR(/i — > 77) is smaller than the SM results due to constraints from the 
conditions of the inflation and the dark matter. These ratio is at most 10 % because our 
model contains only one charged scalar field. 



V. DISCUSSION AND CONCLUSION 



In this Letter, we have not explicitly discussed baryogenesis. It is likely not difficult to 



complement the mechanism for baryogenesis to our model via leptogenesis 47j . In Ref . 



the possibility of the leptogenesis in the Ma model [20J has been studied in details under 
the constraint of current neutrino and dark matter data. By using the typical value for A5 
in our model (A5 ~ 10 -3 -10~ 2 ), the scenario of baryogenesis through the leptogenesis would 
be possible if masses of the right-handed neutrinos are about 10 8 GeV. We will study this 
scenario in more details with its phenomenological consequences elsewhere. 

On the other hand, the possibility of electroweak baryogenesis would also be interest- 
ing Q|. The condition of strong first order phase transition is compatible with rrih = 



126 GeV in the framework of two Higgs doublet models 49( including the inert doublet 



model 



50| | . In such important phenomenologica 



consequence is a large deviation 



in the loop-corrected prediction on the hhh coupling [5JJ], by which the scenario can be 
tested when the hhh coupling is measured at future colliders such as the ILC or the CLIC. 
However, in the inert doublet model including the model we have discussed in this Letter, 
an extension has to be needed in order to get additional CP violating phases, which are 
required for successful baryogenesis. 

We have studied the simple scenario to explain inflation, neutrino masses and dark mat- 
ter simultaneously based on the radiative seesaw model with the Higgs inflation mechanism. 
We find that the parameter region where Z2-odd scalar fields can play a role of the infla- 
ton is compatible with the current data from neutrino experiments and those of the dark 
matter abundance as well as the direct search results. This scenario predicts a specific 
mass spectrum for the scalar fields, which can be measured at the LHC and the ILC with 
y/s = 500 GeV. Our model is a viable example for the TeV scale model for inflation (and 
neutrino with dark matter) which is testable at collider experiments. 
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